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ABSTRACT 

The curvature drift instability has long been considered as a viable mechanism 
for pulsar radio emission. We reconsidered this mechanism by finding an explicit 
solution describing propagation of short-wave electro-magnetic waves in a plasma 
flow along curved magnetic field lines. We show that even though the waves 
could be amplified, the amplification factor remains very close to unity therefore 
this mechanism is unable to generate high brightness temperature emission from 
initial weak fluctuations. 

Subject headings: instabilities - plasmas - (stars:) pulsars: general - radiation 
mechanisms: non-thermal - waves 



1. Introduction 



It is generally believed that the coherent pulsar radio emission is generated by some 
sort of instability in the relativistic electron-positron plasma ejected from the pulsar along 
the open field line tube. One of the widely recognized candidat es is the curvature instabil- 
ity associated with plasma motion along the curved field lines. iGoldreich fc Keeleyl (1l97ll ) 
found an electro-magn etic instab il ity o f a set of charged particles moving on a circular 
ring. On the contrary iBlandfordl (119751 ) demonstrated that maser action is impossible in 
the plasma flow along infinitely strong magnetic field unless points of inflection and tor- 
sion are present in the field g eometry. The controversy was resolved by lAsseo et al.l (119831 ) 
and lLarroche fc Pellatl (119871 ) who showed that the instability develops only if the flow is 
sharply bounded; the necessar y cond itions are unlikely to be met in pulsar magnetospheres. 
Zhelezniakov fc Shaposhnikovl ( 119791 ) noticed that the wave amplification is possible even in 



broad, smooth outflows if the inertial drift of particles is taken into a ccount. This cu r vature - 
drift instability was fu rt her investigated by iKaz 



Lyutikov et al.l ( 11999aU bl): IShapakidze et al.l ( 120031 ): lOsmanov et al.l ( 120091 ). 



jeri et all (11 991): 



Luo fc Melrosd ((19921); 



The curvature- drift instability belongs to the class of resonant instabilities, i.e. it devel- 
ops when and if particles move in phase with the wave so that the power of the electric force 
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applied to the particles, E-v, does not oscillate. The necessary condition for the resonance is 
that the wave phase velocity is less than the speed of light. The transverse wave in a strongly 
magnetized plasma does become subluminal; in the electron-positron plasma moving along 
the magnetic field with the Lorentz factor 7 p , the phase velocity of the wave propagating at 



a small angle to the magnetic filed is (e.g. iKazbegi et al.l (119911 )) 



kc ' 4o;^7p' ^ 



where 



Ane 2 n eB 

uj p =\ ; uj c = — • (2) 

V m mc 

Here e and m are the electron charge and mass, correspondingly, n the pair number density. 
The above simple dispersion law is obtained if the wave frequency in the plasma frame is 
below the cyclotron frequency, u <C ^ c /l P , and if one neglects the energy spread of the 
plasma particles. 

Conventional wisdom is that the electron-positron plasma is generated near the pulsar 
polar cap by a highly relativistic primary particle beam accelerated by rotationally induced 
electric field. Such a two-component flow ideally suits to the curvature drift instability 
because a single component could not experience a significant drift, which requires a large 
Lorentz factor, and at the same time to affect sufficiently the wave dispersion properties. In 
the two-component flow, the dense, relatively low energy secondary plasma makes the wave 
subluminal according to Eq. ([Q) whereas the energetic primary beam has a drift velocity 
sufficient to resonantly excite the wave. 

Because the magnetic field in the pulsar magnetosphere is very large, the factor 5 is 
very small therefore the wave velocity remains very close to the speed of light. Since the 
beam is highly relativistic, the resonance condition, u — k ■ v^,, could be fulfilled only if 
the Wcive prop£i§ cites dbt a very small angle to the beam. However, the radiation propagates 
along nearly straight rays (the refraction index is close to unity) whereas the particle beam 
follows the curved magnetic field lines therefore for any ray, the resonance condition could 
be fulfilled only within a small zone where the ray nearly grazes the field line (see Fig. 1). 

Thus the wave is resonantly amplified only at a short length scale. This is in spite 
of the fact that the size of the unstable region could be large, i.e. one can find unstable 
waves at any point within the region. However, any such a wave rapidly leaves the resonant 
range of the angles when propagates in the curved magnetic field. It is the resonant growth 
length that limits the amplification factor of the wave. The curvature instability could 
be considered as a viable mechanism for pulsar radio emission if the amplification factor is 
exponentially large; only in this case small electro-magnetic fluctuations in the flow could 



Fig. 1. — Geometry of the problem. Radiation propagates along nearly straight rays (dotted 
arrow) across the plasma flowing along the curved magnetic field lines (thick arrows). The 
resonance amplification is possible only if the angle between the wave propagation direction 
and the particle velocity is small. This condition could be fulfilled only within a small region 
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give rise to the observed high brightness temperature emission. We will show in this paper 
that this never happens; the amplification factor in any case remains very close to unity. 

The paper is organized as follows. In the next section, we outline the basic parameters 
of the pulsar magnetosphere. In sect. 3, we present preliminary, qualitative analysis of the 
wave amplification by curvature drift instability. In sect. 4, we derive the wave equations 
and introduce the coordinate system appropriate for analyzing the wave amplification. In 
sect. 5, we find response of the medium, i.e. calculate electric currents excited in the main 
plasma flow and in the high energy beam by the wave. The amplification factor of the wave 
is found in sect. 6. 



2. Plasma flows in pulsar magnetosphere 

Let us first shortly outline the basic parameters of the pulsar plasma. Better conditions 
for the curvature drift instability are achieved in the external part of the pulsar magneto- 
sphere where the magnetic field is not too high. Assuming the dipole field, one finds the 
cyclotron frequency at the distance D from the star as 

u c = 1.8-10 10 ^pHz; (3) 

where 5* is the surface field, d = D/r*, r* the star radius and the usual notation A = 10 X A X 
is used. 

The particle density in the primary beam is of the order of Goldreich- Julian particle 
density, = B/(eP), where P is the pulsar period. Therefore the plasma frequency of the 
beam may be estimated as 

< = = -b 1 - 4 

m P 

Note that both the particle density and the magnetic field in the open field line tube vary 
inversely proportionally to the tube cross section therefore the expression (j3J) remains valid 
in any point. 

The Lorentz factor of the primary beam is very high, ~ 10 6 , therefore they experience 
inertial drift when moving along the curved magnetic field lines. The drift velocity is directed 
along the binormal to the field line; 

FxB 

u d = c—; (5) 



where F = 7f,mc 2 / R c is the centrifugal force, R c the curvature radius of the field line. Within 
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the open field line tube, the curvature radius could be estimated as 

R c = v / ^ = 2-10 9 v / 4Pcm; (6) 
where Rl = P/2tt is the light cylinder radius. Now one gets an estimate 

\u d \/c = 8-10- i lb6 dl- 5 B; 1 \p- 1 / 2 . (7) 

One sees that the drift velocity is very small. Therefore the optimal conditions for the 
curvature-drift instability are achieved in the region of the pulsar light cylinder where the 
drift velocity is maximal. 

The secondary electron-positron plasma is generated with a moderate Lorentz factor 
7 P ~ 10 -T- 100. The plasma density is generally assumed to be large as compared with the 
Goldreich- Julian density. Introducing the multiplicity factor, A, defined as the ratio of the 
plasma to the Goldreich- Julian density, one can write 

u p = Xu b = — — . (8) 

According to the ayailable polar cap models, A varies from a few to a few thousands 
(IHibschman fc Arond 120011 ). Obser yations of puls ar wind nebulae give evidence for much 
larger multiplicities, A > few x 10 5 (jde Jagerl 120071 ). 

Now one can estimate the parameter S, which defines, according to Eq. (OQ), the disper- 
sion properties of the electro-magnetic waves in pulsar magnetospheres: 

5=1.7-10~ 6 -M-. (9) 

One sees that the phase velocity of the waves is very close to the speed of light at any 
reasonable pulsar conditions. 



3. Preliminary analysis of the problem 

We are going to study the resonant amplification of radio waves when they propagate 
across the pulsar plasma. As it was already discussed in the Introduction, the wave could 
be resonantly amplified only if it propagates close to the direction of the particle beam. The 
particles move along the curved magnetic field lines and also could experience a weak inertial 
drift along the binormal to the field line. Therefore the amplifying wave should be directed 
at a small angle to the magnetic field and polarized predominantly along the direction of the 
drift velocity. Since the radiation propagates along nearly straight rays whereas the beam 
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particles follow the curved magnetic field lines, the amplification could occur only in a small 
zone (Fig. 1). The magnetic field lines within this zone could be considered as arches of 
concentric circles along which both the plasma and the beam circulate. This configuration 
approximates locally real plasma motion along curved magnetic field lines in the pulsar 
magnetosphere. Therefore one can conveniently use the cylindrical coordinates (r, 0, z) with 
the direction along the background magnetic field. 

In this configuration, one can expand the wave in cylindrical modes, 

E(r, t) = E W;S)fcz (r) exp(-iut + ik z z + is<f>); (10) 

thus reducing the problem to an ordinary differential equation for E WjSj ^(r). The wavelength 
is very small as compared with the characteristic inhomogeneity scale, 

cor/c = 4- IOVghzV^; (11) 

where /gh z is the wave frequency in GHz; we used the estimate (jUJ) for the curvature radius. 
Therefore one can try to use the WKB approximation. 

In this case, the solution should take a form of locally plane waves satisfying the dis- 
persion equation ([[]), i.e. there should be 

E(r,t) oc A out e^ krdr + A ia e~ i f krdr ; (12) 

where the radial wave vector, k r , should be found from the dispersion equation (1) as 

cu^ s 2 

k " = c 2 (l - 5f (13) 

The first term in Eq. (fT2l) represents rays propagating away from the axis. At any point, 
the direction of such a ray is along the wave vector (k r , = s/r, k z ). Taking into account 
Eq. ( I13p . one finds the impact parameter of the ray, 

= ; C ! (14) 



We are interested in rays propagating together with the beam therefore k z /u ~ < 1. 
Then one can expand 

cs f 1 c 2 kf 



<'o = -(l-<5 + ^l- (15) 



The second term in ([T2"j) represents rays propagating towards the axis with the impact 
parameter (IToT) . Therefore the cylindrical wave could be considered as a set of rays with 
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the same impact parameter (Fig. 2a). These rays converge to the axis till the minimal 
distance r and then go to the infinity. Two rays, ingoing and outgoing, pass through any 
point at r > r . At r « r , a caustic is formed where the WKB approximation is violated. 
The structure of the solution should resemble that of the Airy function (Fig. 2b): the 
maximum near tq, oscillations at r > r$ and the exponential decay (the evanescent wave) at 
r < r$. The width of the caustic is very small, Ar = r^ 3 (c/uj) 2 ^ 3 <C tq. The oscillating part 
approaches the WKB solution (fT2"l) at r — r ^> Ar. 

In the r — plane, the lines of constant phase rotat e with a constan t angula r velocity u/s. 



At r > ro, they form unwinding spirals (see Fig.2 in iLyutikov et al.l (11999 bl ) ) approaching 
eventually the lines J k r dr + s<fi/r — const. Normals to these lines form rays plotted in Fig. 
2a. The azimuthal component of the wave phase velocity f p h,</> = w/A^ = ruj 7 ' s var ies with 
the radius. Therefore the resonant interaction with particles could occur only in a narrow 
layer where the particle azimuthal velocity is close to ru/cs. 

More exactly, the resonant condition is u = k • v. The particle azimuthal velocity is 
Vfj, = ^Jc 2 —u 2 d — 1/(27^) ~ c [l ~~ M d/(2c 2 )] therefore the resonance condition is written as 

cu = k z u d +(cs/r){l-u 2 /(2c 2 )]. (16) 

Resolving this equation with respect to r, one finds the radius of the resonance layer 

(17) 

The condition that the wave reaches the resonance radius, r < r res , can now be written as 



cs 



1 | u d ( ck z u d 
c \ u 2c 



f lLvutikov et al.lll999ah 



2 \ UJ c 

The right-hand side of this expression goes to zero if the ray is inclined to the magnetic 
field line plane at the same angle as the particle beam (more accurate calculation shows that 
the right-hand side goes not to zero but to % 2 , which could be neglected) therefore at first 
glance, this condition does not impose any restrictions on 5. However, we will see that the 
amplification also ceases in this case therefore some discrepancy between the directions of the 
ray and the beam is necessary. The reason is that the resonance interaction occurs via the 
component of the wave electric field parallel to the particle velocity, E • v. The transverse 
wave should be inclined to the beam in order for this component be non-zero. Then Eq. (fl~8l) 
implies that the curvature drift mechanism could operate only if the magnetospheric plasma 
is dense enough. 

Since both 5 and u d /c are very small at any reasonable pulsar conditions, the ray should 
be inclined to the plane of the magnetic field line by a small angle. In this case, the resonance 
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and the caustic radii are close to each other. Therefore in order to find the amplification 
factor of the wave, we could solve the wave equations in a relatively narrow range of radii 
containing both r and r res . At the external boundary of this region, the solution should 
be matched with the WKB solution ( lT2"j) . At the internal boundary, the solution should be 
matched with the evanescent wave. Then the amplification factor is found as the ratio of 
amplitudes of the outgoing and ingoing waves, a = \Ai n \/\A ont \. With such a setup, we have 
in fact to solve the reflection problem for the radial wave E w jS) fc z (r). Therefore the curvature 
instability may be thought of as a part ial case of the ge neral effect, namely, overreflection of 
waves from a rotating dielectric body ( jZerdovichl 119721 ) . 



4. Wave equations 

The propagation of the wave is described by Maxwell's equations 

V-E = 4vrp; V x B = 4vrj + § ; (19) 
V-E = 0; VxE = -f; (20) 

where currents and charges are excited in the plasma by the wave. From here on we take 
the speed of light to be unity. One can conveniently study the curvature instability by 
expanding the wave in cylindrical modes ( flOl) . Eliminating B, one writes the equations for 
the amplitude of a cylindrical wave as 

d 2 Es 1 dEs ( n , 9 1 \ _ s d E r k z s _ . „. 

-TT + --T+ W-k\-- )E 4> -i- — ^ + ^E z = -Ancutjf, 21 
ar r dr \ r J r dr r r 

d 2 E z ldE z f n s 2 \ ^ k z d ^ k z s ^ . . 

-r^ + --r^+ (cu 2 --)E z -i^—rE r + ^E^ = -4iruij z] 22 
ar z r dr \ r A J r dr r 

s 2 \ d dE 

tu 2 — k 2 J E r — is—rEj, — ik z —^- = —4iruij r ; (23) 

r 2 J dr dr 

We are going to solve these equations close to the caustic zone where the wave propagates 
nearly along the direction (which coincides with the direction of the background magnetic 
field). We will see that the wave remains nearly transverse therefore one can conveniently 
introduce the variables 

E+ = JL E + t± Ez] E _ = ± Ez -^E (24) 
ur uj ur uj 

such that in the region of interest, E + becomes the longitudinal component of the field 
whereas E_ and E r the transverse components. Combining Eqs. ( 1211) and (|2"2"|) . one gets the 
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equation for EL 
d 2 E_ 



1 I i 
+ - 1 + 



o?r 2 r \ u; 2 r 2 



LO 



1 + 



w 2 r 2 



E_ 



/c r s dE 4 



2sk„ 



where 



3- 



Lo 2 r 2 dr Lor 



s . k z . 

Jz 3<j>- 
Lor lo 



2 



-4niLoj. 



We complement this equation by Eq. fl23|) for E r and by the Gauss law 

ldrE r 



r dr 



+ iuE + = 47cp. 



(25) 
(26) 

(27) 



Sin ce ur a s > 1, we can solve these equations in the short wave approximation (e.g. 



Nayfehl (119731 )). The resonance occurs close to the caustic zone therefore one can conveniently 
use the caustic coordinates. The plasma refraction index is very close to unity therefore the 
characteristic width of the caustic zone is the same as in the vacuum case, Ar ~ r/s 2//3 ; 
the reflection point is just shifted towards a smaller radius. Therefore let us define the new 
radial dimensionless coordinate, 



r [l + 



x 

^2/3 



(28) 



where r is the WKB impact parameter (IT5|) . which is in fact the caustic radius. Transforming 
to the new variable and retaining only the leading order in s -1 / 3 « (wro) -1 / 3 terms, we can 
reduce the equations f l25|) . fl23l) and (1271) to the form 



s 2 / 3 5)E_ 



E"_ + 2(x 

is 1/3 E' + -2(x-s 2/3 5)E r = 

E' r + is 1/3 E + = 
where prime denotes differentiation with respect to x. 



4ns 2 / 3 
4ns 2 ' 3 
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47TS 1 / 3 

/ 

LO 



-I3r 



(29) 
(30) 
(31) 



In order to write the equations in the closed form, we have to find the charge and current 
densities excited by the wave in the plasma and in the beam. Here we consider the simplest 
case when both the plasma and the beam are cold so that they are described only by their 
densities and velocities. Then the current and the charge densities are presented as 



P 



E 



q a n c 



E 



(32) 
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where the summation is over the particle species (electrons and positrons in the plasma and, 
say, positrons in the beam). Now we have to find perturbation of the particle motion in the 
electro-magnetic field of the wave. 



5. The response of the medium 

In this section, we find the current and charge density induced in the medium by the 
wave. With this purpose, we have to solve the equations of particle motion in the field of 
the wave: 

-J— + (v V) 7 v = — (E + v x B + v x B). (33) 
at m 

Here q = ±e is the particle charge; the background field, B , is assumed to have only <fi 
component. The oscillatory motion of the particles in the electro- magnetic field of the wave 
implies perturbation of the particle density according to the continuity equation 

— + (V • nv) = 0. (34) 

Without the wave, the particles move in the <fi direction along the background magnetic 
field and also experience the inertial drift in the z direction. In order to find the response 
of the medium to the wave, one has to linearize Eqs. ( 133|) and (|34|) with respect to small 
perturbations of the velocity, v = + and of the density n = n^ ' + n^'. The response 
is expressed via the solution to the linearized equations as 



P 



E rf; j = E ^M 0) + -L 0) v«). (35) 



Even though the procedure is straightforward, it leads to rather cumbersome expressions 
therefore we consider response of the plasma and of the beam separately exploiting from the 
beginning appropriate approximations. 



5.1. Response of the beam 

The Lorentz factor of the beam is very large, 7& ~ 10 6 , therefore when moving along 
the curved magnetic field line, it experiences the inertial drift in the z direction 

^,beam = U d = ~— ~ • (36) 

The unperturbed velocity in the direction may be presented as 



«$-» = VWl-ri2- 76 - 2 . (37) 



- 11 - 



In this subsection, only velocities and density of the beam appear therefore we drop the 
index "beam" in v*^ and rS l \ Linearizing the continuity equation f )34|) one gets 

where 

Q — u — V- — k z u d + iO; (39) 
r 

and we take into account that according to the Landau prescription, a small imaginary part 
should be added to the frequency. 

The linearized equations of motion are written as 

^(D + %(i) + ^W = J-(VB,-E r -u i Bfr (40) 



lb r r m^ h 

Vu d E z - \u 2 d + ^) E (f> + u d B r 



r 7717& 



(41) 



iilvP-^-v^ = -—[(l-u 2 d )E z + VB r -u d VE^ (42) 

lb ™lb 

The terms with V/r in the left-hand side of Eqs. ( l4"0l) and ( T4T]) arise from the centrifugal 
and Coriolis forces, correspondingly. These terms were neglected in all previous works on 
the curvature instability. We will see that they are important near the resonance where Q 
goes to zero. 

The density of the beam is extremely small therefore contribution of the beam to the 
current and charge densities of the medium is significant only near the resonance layer where 
the wave moves together with the beam so that where Q goes to zero. The resonant radius 
is defined from the condition Q = 0, which is nothing more than the condition 00 = k ■ V5 
already discussed in Sect. 3. Therefore we can use Eq. (|T7|) for r rcs . Expressing r res via 
the radial coordinate x defined by Eq. (j2"81) and making use of Eq. ffTBT) for r , one gets the 
position of the resonance layer in the dimensionless variable 

Zres = 5* 2/3 - \s 2 '\l - K)W d . (43) 

Taking into account that the ray is directed nearly along the particle beam, one can expect 
k z jbj ~ u d therefore we introduced the parameter 

k = (44) 
u d cu 

The function Q can now be written via the dimensionless variable as 

n = - x vcs + i0). (45) 
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When solving Eqs. ( HUH21) for in and we retain only resonant terms (those 
proportional either to Q~ l or d/dr). Expressing B via E with the aid of Maxwell's equations, 
one finds 



«W = -i 



1 



eu d 



1 n ^ 

— — K{1 — K)U 



ui d 



m^V 2 1 Vt 



— n(l — n)u d 



-uV^ {VE " +U ' lE): (1C) 
+ } {VE+ + u d E z ) . (47) 



The component Vr has no resonant terms however, it should be taken into account in the 
expression (138]) for because being differentiated, it could contribute to the resonance 
interaction provided the field amplitudes are singular. Dropping the terms proportional to 
Q one finds 

v V = _ k) {ve + ^ (4g) 

muj c 

One sees that the resonance interaction is determined by the projection of the wave elec- 
tric field onto the direction of the beam. This projection is expressed via the "longitudinal" 
and "transverse" components of the electric field (124"]) as 

VE, + u d E z =( S ^ + ^)e + +(^-^-)e^E + + (1- n)u d E_. (49) 

Substitution of Eqs. (146]) . (|4"Tj) and f )48|) into Eq. fl38|) yields the beam contribution to 
the charge density in the form 

4™ P bcam = Anient = - ^^^2 f 1 + ( x ~ x ^4~) i E + + ^ ~ K ) U * EJ \ • ( 50 ) 

Here we made transformation to the dimensionless coordinate x and also dropped the term 
1/7^ as compared with u\ in the square brackets in Eqs. ( 1461) and ( j4"Tl) . One sees that the 
density perturbation is of the second order in Vt~ l oc (x — x rcs ) _1 whereas the velocities are 
only of the first order. Therefore the beam contribution to the current could be written as 

;beam _ beam (0) 

J — P V beam> \ 01 ) 

which yields 

^-bearn = ^beam (Ky _ J_ U A ~ (J _ K )u d p hea ™. (52) 



oj tur 

Note that if we were neglected the V/r terms (arising from the Coriolis and centrifugal 
forces) in the linearized equations of motion (14*0]) and fj4"T|) . we would come to the expression 
similar to ( 150]) but with the coefficient proportional to u d instead of u\. Therefore taking 
into account these terms is crucially important. 
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5.2. Response of the plasma 

The Lorentz factor of the plasma particles, j p — (1 — f 2 )" 1 / 2 , is not very large so that 
one can neglect the inertial forces. Then the unperturbed velocity is in the <fi direction, 
vi°^ = 0; v? = v p . In this subsection, we drop the index "plasma" in v^ 1 ) and The 
linearized equations of motion are now written as 

^-vK ,+ ^" = < 53 > 

iL-^)vf = (54) 
V r / v 

z(u;-^U«-%« = -^-(E z + v p B r ). (55) 

Since the plasma is electrically neutral in the sense that the densities and unperturbed 
velocities of electrons and positrons are equal, the plasma response is determined only by 
terms with odd powers of q in the perturbed velocity These terms are easily found as 



(1) _ ej p (u-v p s/r 

U T . V 



muolixi 



[u — v p s/r)E r — iv p - 



+ {terms oc e 2 }; (56) 



dr 

vf = T\ E * ^ 

v m^p[oj — VpS/r) 

= -i eip ( U ~ *W r ) [(^ _ VpS / r )E z + v p k x Ej\ + {terms oc e 2 }. (58) 

TflU) c Ul 

Here we used the strong field approximation, 

w c > jp(uj - svp/r) « w/27p. (59) 

The charge and current densities are written as 

pP ia Sma = 2en (i) = Zerip f s m + ^ _ zdr^A 
— f p s jr \r 9 r dr J 

it = «pP + 2en4 1} ; * = 2en p T;< 1 >; j r plasma = 2en p t;< 1 >. (61) 

Here n p is the pair number density (the number density of electrons or positrons); the factor 
2 takes into account that contributions of electrons and positrons are equal. Note that 
the plasma current and charge densities are small as 5 ~ (u p /u c ) 2 (see Eq. (JTJ). Then 
the "longitudinal" component of the electric field, E + , is also small as 5 (or as s -1 / 3 , see 
Eq.f l3ip ) therefore one can neglect the contribution of this component to the plasma current 
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and charge densities. Then one gets 



plasma <-> J ^p^P 



4mf_~ = 2 



,2, 



U! 3 U! 2 L 



US VpS 



1 



u 2 p k 2 z \ E 2 U 2 p k 2 z d 2 E_ | u$k 2 dE r 
^(cu — Vps/r) J u 3 u 2 dr 2 ^ r 

Here we took into account that dE^/dr ^> E_/r. 

These expressions could be significantly simplified taking into account that we are inter- 
ested in the region close to the resonance layer defined by Eq. (17) so that we can substitute 
r res instead of r. Then one can write, e.g., 

sv p { 1 — KU 2 d \ UJ 



In the last equality, we take into account that 7 P |wd| "C 1. Making use of this approximation 
and taking also into account the condition (I59I) . one finds 

A • -Plasma 2 f U 8 7p^ \ p , U P k * f ; U p/ 2 ^ P /A 

The term with E' r can be neglected because one can show with the aid of Eq. ( 13T|) that this 
term is small as compared with the term with E_ (small as E + /E_ or as uE_/p ~ 5). The 
term with E 1 " could also be neglected because when substituting the expression for jP lasma 
into the right-hand side of Eq. ( |29|) . one sees that this term is small as compared with the 
corresponding term in the left-hand side. Assuming also that k z /u ~ Ud <C u;/(4a; c 7 2 ) one 
finally finds 

47r2J plasma = = 2 6u)E-. (65) 

As a consistency check, let us show that this expression for the plasma response leads to the 
the dispersion law JT]). 

With this purpose, let us substitute this expression into the wave equation ( I2"§j) and 
return to the original variable r. This yields 

r d2E ~ , o_2/3 r - r 



s 4/3 ^ r 2 rQ 

In the WKB approximation, £L = exp(i J k r dr), one gets 



+ 2s 2 / 3 - = 0. (66) 



s 2 

A; 2 = 2— (r - r ). (67) 



On the other hand, the dispersion equation (P) is written in the form ffl3l) . which is reduced 
to flBT|) at r — r r . 
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6. Amplification factor 

6.1. The wave equation in the closed form 

With the results of the previous section, we can write the wave equations in the closed 
form. Substituting Eqs. f )50|) . f )52|) and fl65|) into the equation for the "transverse" component 
of the electric field f )29|) one gets 

E"_ + 2xE_ = a (l + {x- x Ies )-f) ( £L + E+ ) ; (68) 

(x-x ICS + tOy V dxj V {l-K)u d J 

where 

Taking into account that s ~ cor, one can present a, with the aid of Eqs. (BJ, (EI) and (136]) . 

as 

a = 2(D/R L ) l/2 u 5 d K(l - k) 3 . (70) 

Taking into account that D < R^, n ~ 1 and the drift velocity is always small (recall that 
here Ud is measured in units of c), one sees that a is a very small quantity. This already 
ensures that the amplification factor is small but let us calculate it explicitly. 

Eq. (168]) should be generally complemented by Eqs. (130]) and (131]) . However, Eq. (1681) 
is reduced to a closed equation for E_ provided 

E + < (1 -K)u d E„. (71) 

The conditions for such an approximation could be found as follows. 

One can estimate E + from Eq. (1311) . Outside the resonance region, one can write 
E + ~ s-^Er because the right-hand side of this equation is of the order of SE_ s 1 ^ 3 E_. 
We are interested in the wave polarized predominantly in z direction therefore E r ^ E-. 
Then the condition flTTl) is fulfilled if 



(1 - K)s l l 3 u d > 1, (72) 

i.e. if the drift velocity is not too small. Taking into account that the effect under con- 
sideration assumes non- negligible u d so that the wave amplification is suppressed when u d 
goes to zero, we assume that the condition (1721 is fulfilled, which justifies the inequality (17T]) 
outside the resonance region. It is shown in Appendix that if the inequality (171]) is fulfilled 
outside the resonance region, it remains valid also within the resonance region. Therefore 
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the condition f J72|) justifies neglect of E + in Eq. fl68|) . Now we come to the closed equation 
for E_ 



E"_ + 2xE_ 



a 



(x - x res + ioy 



E_ + (x 



dx 



(73) 



Inasmuch as a < C 1, thi s equa tion could be solved by the method of matching asymptotic 



expansions (e.g. 



Navfehl (Il973^ . 



6.2. Solutions outside and inside the resonance zone 



At \x — x res | ^> y/a, one can neglect the term the right-hand side of Eq. ( 1731 ; then the 
solutions are expressed via the Airy functions. Specifically to the left of the resonance, at 
x res — x ^> y/a, we have to choose the evanescent wave solution 

^nleft 



= Ai(-2 1/3 0- (74) 
3> -y/a, the external solution should be chosen in 



To the right of the resonance, at x 
the general form 

E" ght = A l M(-2 1 ^) + A 2 Bi(-2 1/3 £); (75) 

the amplitudes A\ and A\ should be found by matching the solution (1751) to the solution 
fl?4"j) via the resonance layer. In the wave zone, x ^> 1, the function (1751) is reduced to a 
superposition of an ingoing and outgoing waves 

ciright _ 1 

2 7 /6 x l/2 

The amplification factor is defined as the ratio of the amplitudes of the outgoing and ingoing 

waves 

Ai + iA 2 



(A, + iA 2 )A^ /2 ^ ] + (A, - zA 2 )e-^(^) 3/2 -f ] 



(76) 



A 1 - iA 2 

We will see that \A 2 \ <C \A\\ (at a = 0, the full solution is described be Eq. ( JT4"j) . which 
yields A\ — 1; A 2 — 0; therefore at a < 1, we obtain A 2 <C Ai). Then one can write 



1 + 2%A 2 



(n 



so that the amplification factor is determined only by the imaginary part of the amplitude 
A 2 . Now let us proceed to matching the solutions _E" ght and E l ^ ft via the resonance zone. 

First we have to find the "inner" solution, i.e. that in the resonant zone \x — x res \ <C 1. 
Introducing the "inner" variable z = (x — x r 

d 2 E 



a, one writes Eq. fl73|) as 



dz 2 



+ 2a(y/az + x res )E. 



[z + iOf 



E_ + z- 



dE_ 

dz 



(79) 
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One sees that one can neglect the second term in the left-hand side; then the solution is 
easily found in the power-law form. Returning to the original variable, x, one can write the 
solution in the resonant zone as 

E_ = Q(x - x res + i0)~ a + P(x - x rcs + i0) 1+2a ; (80) 

where Q and P are constants. Now we have to match the solutions (|74p . (|75p and f[8"Uj) . 
This could be done because the "outer" solutions (174"|) and ([75]) are valid at x rcs — x ^> y/a 
and x — x res 3> \/a, correspondingly, therefore at small a, the domains of validity of these 
solutions are overlapped with the domain of validity of the "inner" solution, \x — x rcs \ <C 1. 



6.3. Matching the solutions 

We first find the constants Q and P by matching the "inner" solution ( IBUl) with the 
solution (|74p . In the region 

\fa <C x res - i C 1 (81) 



the solution ( 1801) is still valid. Having found Q and P one can find amplitudes A\ and A2 
by matching the solution fjHUj) with f }?4"j) in the region 

y/a < x - x rcs < 1. (82) 

We will see that it is not necessary to perform the full matching procedure if we are interested 
only in the reflection coefficient, which is determined only by the imaginary part of the 
amplitudes. 

In the region ( 1HTT) the "inner" solution ( 180]) is reduced to 

E_ = Q (1 — a In \x — x ms \ — ma) + P(x — a; rcs ) (1 + 2a In \x — x res | + 2ma) . (83) 

The "external" solution is found from Eq. (173"j) as a slightly perturbed evanescent wave. 
Namely the right-hand side of this equation is small in the region (I8"T|) therefore one can 
substitute there the evanescent wave (jT4|) thus coming to an inhomogeneous Airy equation. 
The solution to this equation, which goes to (174]) at x — > — 00, is written as 



E_ = E Mt + a 



x G(x,xi) 

00 (^1 ^res) 



1 + [Xi - x res ) — 
axi 



E^ix^dxt; (84) 



G(x, Xl ) = 2 1 / 3 tt [Ai(-2 1 / 3 x 1 )Bi(-2 1 / 3 x) - Bi(-2 1 / 3 x 1 )Ai(-2 1 / 3 x)] . (85) 
At x — > x res the first term in this solution is reduced to the linear function 

E Mt (x -> a; rcs ) = Ai(-2 1 / 3 x rcs ) - 2 1 / 3 Ai'(-2 1 / 3 a; rcs )(x - a; res ); (86) 
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so that matching with the "inner" solution fl83l) to within the zeroth order in a yields 

Q = Ai(-2 1 / 3 x res ) + 0(a); P = -2 1 / 3 Ai'(-2 1 / 3 x res ) + 0(a). (87) 

As x — > x res , the second term in (1841 gives the terms of the order of a In | x — x res \ and a 
therefore the solution (1841) could be smoothly matched with the "inner" solution (183]) . 

We are interested only in the imaginary parts of the coefficients Q and P, which are 
sufficient in order to find the imaginary part of the amplitude A 2 . Since the solution ( 184"]) is 
real, it is matched with the solution (183]) if = na^RQ and QP = —2-Ka^tP to within the 
first order in a. Making use of Eq. (1ST]) , one gets 

%Q = -7raAi(-2 1/3 x res ); = -2 4/3 7raAi'(-2 1/3 x res ). (88) 



In the region (18 2 p the resonance solution (18 Op is written as 



£L = Q [1 — a \n(x — x res )} + P(x — x TCS ) [1 + 2a \n(x — x res )] ; 
whereas the "external" solution is presented as 



E_ = E 



right 



a 



At x — > x res , one can substitute E_ by a linear function 



G(x,x 1 ) 

%1 %res) 
right 



/ \ d 

1 + (xi -x res ) — 
ax\ 



E r y ht (x 1 )dxi, 



E 



right , 



x) = A!Ai(-2 1/3 x res ) + A 2 Bi(-2 1/3 
-2 1 / 3 [A 1 Ai'(-2 1 / 3 x res ) + A 2 Bi'(-2 1 / 3 a; rcs ] (x - 



(89) 



(90) 



(91) 



The second term in Eq. (190]) has the order of aln(l/a) in the region ( 182]) because the integral 
diverges only logarithmically. Therefore matching of the solution (I90p with (18"9~]) in the zeroth 
in a approximation yields, with account of Eq. (18T]) . 



A 1 = 1 + 0(a) 



A 2 = 0(a) 



(92) 



In this case, the general expression for the amplification factor (177]) is reduced to Eq. (178]) . 
Therefore it suffices to find only the imaginary part of A 2 . 

With this purpose, we can only match the imaginary parts of (1891 and ( 190]) at x — > x res . 
With account of Eq. (188]) . one gets from Eq. 



%E_ 



-na [Ai(-2 1 /3 Xrcs) + 2 4 / 3 Ai / (-2 1 /3 Xrcs )( x _ Xres) ] + (a 2 ). 



(93) 
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Eq. fl90l) yields at x — > x vcs , with account of Eqs. f l9T]) and f )92|) . 

3£L = SA!Ai(-2 1/3 x r c S ) + 5M 2 Bi(-2 1/3 x rcs ) (94) 
-2 1 / 3 [SA 1 Ai , (-2 1 / 3 x res ) + SM 2 Bi'(-2 1/3 x rcs ] (x - x rcs ) + 0(« 2 ). 

Comparing the coefficients at (x — x res )° and (x — x rcs ) in the two last equations one gets 

ZAi = -«vr 2 [Ai(-2 1 /3 Xrcs )Bi'(-2 1 /3 Xrcs ) + 2 4 / 3 Ai / (-2 1 / 3 x 1 . cs )Bi(-2 1 / 3 x res )]; (95) 
%A 2 = (l + 2 4 / 3 )avr 2 Ai(-2 1 / 3 x res )Ai , (-2 1 / 3 x res ). (96) 

Here we used the formula for Wronskian, AiBi' — Ai'Bi = 1/ir. Now Eq. (178"j) yields finally 

a = 1 + 2 (1 + 2 4/3 ) 7r 2 aAi(-2 1/3 x rcs )Ai / (-2 1/3 x rcs ). (97) 

One sees that the wave is amplified, a > 1, provided Ai(— 2 1/,3 x res )Ai'(— 2 1//3 x res ) > 0, which 
happens if x res falls within the intervals 



1 < 2 1/3 x TCS < 2.4; 3.2 < 2 1/3 x res < 4; 



3 



l + 4Jfe) 



7T 



2/3 



<C 2 ^ x rPS <C 



3 



-1 2/3 



;(3 + 4fc)7T 

1 

(98) 

where k is a large integer number. 

It follows immediately from Eq. fj4*3"]) that the necessary condition for the amplification, 
x res > 2 -1 / 3 , is achieved only if 5, i.e. the deviation of the wave velocity from c, is not too 
small: 

Taking into account the condition (!72|) . one sees that the condition (|99|) roughly coincides 
with the condition (TIB"]) found in the WKB approximation. This condition implies that the 
curvature drift instability is possible only if the pulsar plasma density if high enough. 



In any case one sees from Eq. fl97|) that even the instability condition is fulfilled, the 
amplification factor f )97|) is very close to unity because it follows from eq. f JTOj) that a 
remains small at any reasonable conditions. Therefore the curvature drift instability could 
not amplify small fluctuations up to high brightness temperatures characteristic for pulsar 
radio emission. 



7. Conclusions 



In this paper, we rederived the curvature- drift instability, which has long been considered 
as a viable mechanism for pulsar radio emission. Our approach differs from that adopted 
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in the previous works on the topic. Namely we did not look for the growth rate of the 
instability but instead explicitly calculated propagation of the electromagnetic wave through 
the plasma moving along the curved magnetic field lines. With such an approach, we found 
the amplification factor of the waves. In the standard approach, the amplification factor 
depends on the growth rate and also on the size of the resonance zone. In the case under 
consideration both the growth rate and the resonance zone are determined by the curvature 
radius therefore there is no sense to calculate both these quantities separately. Moreover, 
in the previous works the growth rate was found by the dielectric tensor method, which 
implies that the waves are considered as locally plane. Then one could define the dielectric 
tensor of the medium such that the growth rate is an imaginary part o f an appropriate 
eigenvalue. Even in the most advanced work by iLyutikov et al.l (Il999bl laf). where it was 
explicitly demonstrated that the resonance occurs in the vicinity of the caustic zone so that 
the wave should be described by the Airy type functions, the amplification was calculated 
only at the WKB condition x res 3> 1. Our approach is free from these limitation and 
moreover provides the amplification factor straightforwardly without intermediate steps. 

The advantage of explicit solutions over the diele ctric tensor approach was clearly 
demonstrated by the long-lasting discussion initiated by iBeskin et al.l (119871 . Il988al ) works. 
These authors derived the dielectric tensor of a weakly inhomogeneous plasma in the in- 
finitely strong magnetic field and found unstable modes. The co r rectness of the obtaine d 
dielectri c tensor has being ext e nsively debated, see Nambu ( 1989 ); Machabeli (jl991 . 1995 ); 
Istominl Jl994j ): iNambul Jl996h: iBornatici &: Kravtsov! (120001 ) . However, the claim of Beskin 
et al was disproved by lLarroche &: Pellatl (119871 ) who presented an explicit solution to the 
wave equations and demonstrated that t he wave ampl i tude does not grow unless t he flo w 
is sharply bounded (see also discussion in IBeskin et al.l (Il988bl ); lLarroche fc Pellatl (119881 )). 
The specific property of the obtained solution is that it becomes singular in a narrow reso- 
nance layer. The dielectric tensor method assumes implicitly that the global solution could 
be described by a locally plane wave with smoothly varying parameters. In shear flows, there 
are typically no such solutions satisfying reasonable boundary conditions; the wave becomes 
singular in the critical lay e r whe re the phase velocity becomes equal to the flow velocity (e.g. 
Stepanyants fc Fabrikantl (119891 )). Therefore the dielectric tensor method does not provide 
unambiguous conclusions on the stability of shear flows; one has to find explicit solutions as 
it was done in the present paper. 

One more novelty of the present research is taking into account of the inertial forces 
(Coriolis and centrifugal) when considering interaction of the waves with resonant particles. 
These terms were ignored in all previous works on the curvature drift instability. We have 
shown that these forces significantly suppress the growth rate. 
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The net result of our study is that at no condition the amplification factor could become 
large therefore the curvature drift instability should be excluded from the list of potential 
mechanisms for pulsar radio emission. 

This work was supported by German-Israeli Foundation for Scientific Research and 
Development under the grant 1-804-218.7/2003 and by the Israeli Science Foundation under 
the grant 737/07. 



A. Appendix. General solution to the wave equations in the resonance region. 

Within the resonance region, \x — x res \ <C 1, one has to retain only the resonant terms 
so that the equations (12U|) . (1311 and (}3"Tj) could be written, with account of (1501 and (152"j) . in 
the form 

E ~ (x - x TCS + iO) 2 ( 1 + ^ Xrcs ^ dx) + (1 - K)Ud) ' 
is 1/3 E' + = 2x Tes E r ; (A2) 

a-w* - i {x _ + ffl) , (i + (» - (g- + ■ (A3) 



Eliminating E r from the last two equations and making use of Eq. (f43|) . one gets 

(1 - .Ka / d\ f E _ + E+\ 

(x - x rcs + «0) z V rfxy V (i - «)W 

thus reducing the system to two equations for .EL and EL. 

It follows immediately from Eqs. (jA"Ij) and (TAij) that + (1 - K)u d E"_ = 0, which 
implies 

J E+ + (l-K)u dJ E_ = C 1 + C 2 (a;-x res ); (A5) 

where and C*2 are constants. Now one can find the fields in the resonance region in the 
form 

E_ = Q + P(x- x Ies ) + a[S{x ~ + 71 ln(x - x res + iO); (A6) 

(1 - KjUrf 

£+ = -T - (1 - K)u d Q + ^ [5 - 2(1 - K)u d P] (x - x res ) 

-a{S(x - x rcs ) + T] ln(x - x res + iO); (A7) 

where P, Q, S and T are constants. These constants could be found by matching, at 
\ x — x tcs\ ~ 1; with the solution outside of the resonance region. 
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It has been shown in sect. 4 that at the condition (j72p . the "longitudinal" component, 
satisfies the inequality ( 171]) outside the resonance region. Inspecting Eqs. (IA6[) and ( 1A7|) 
one sees that this inequality is fulfilled at \x — x TCS \ ~ 1 provided 

|T + (l-«)u d Q|<(l-«)u d Q; (A8) 

\3 - 2(1 - K)u d P\ < (1 - K)w d -P; (A9) 

which yields 

T = —(1 — K)UdQ] S = 2(1 -k)u d P. (A10) 

One immediately sees that in this case the inequality (ITT]) is satisfied also within the resonance 
region with the exception of the negligibly small region \x— x res \ < exp(— 1/a). Therefore one 
can neglect E + in the resonance term (that in the right-hand side of Eq. (I68p ) thus coming 
to the closed equation (175]) . This justifies use of Eq. ( 175]) in order to find the reflection 
coefficient. Note that upon substituting Eq. (lAlOj) into the general resonance solution ( IA6j) 
one comes to Eq. (189]) obtained as the solution to Eq. f]73]) . 
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Fig. 2. — The structure of the cylindrical wave, a) In the WKB approximation, the wave 
is comprised from a rays (thin arrows) coming from all directions with the same impact 
parameter, r . At r = r they form a caustic (thick circle), b) The radial structure of the 
wave in the caustic zone. The vertical stripe shows the position of the resonance layer. 



